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ABSTRACT 


The simulation experiment described herein addresses the effects of structural flexibility 
on the dynamic characteristics of a generic family of aircraft. The simulation was performed 
using the NASA Langley Visual/Motion Simulator facility. The vehicle models were obtained 
as part of this research project. The simulation results include complete response data and 
subjective pilot ratings and comments and so allow a variety of analyses. The subjective 
ratings and analysis of the time histories indicate that increased flexibility can lead to 
increased tracking errors, degraded handling qualities, and changes in the frequency content of 
the pilot inputs. These results, furthermore, are significantly affected by the visual cues 
available to the pilot. 


CHAPTER 1 
INTRODUCTION 


With the advent of larger aircraft and the ever-present desire for lighter structures for 
improved fuel economy, there is increased potential for significant contributions from 
aeroelastic effects in the dynamic response of flight vehicles of the future to pilot inputs and/or 
atmospheric turbulence. Earlier results^ 1,2,3 ! from the present research program indicate how 
aeroelastic effects can be significant and that, as stiffness is reduced, these effects become 
magnified. Not only can the dynamic response due to the aeroelastic modes themselves 
become significant, but coupling between the aeroelastic and the rigid-body degrees of 
freedom can noticeably alter the “rigid-body” dynamics of the vehicle. ^ 2 ’ 3 ^ 

One of the few generic experiments in the literature studying aeroelastic effects on pilot- 
vehicle performance was conducted by YenJ 4 ^ Yen’s experiment, using a very low bandwidth 
longitudinal task, exposed the potential for significant degradation in performance due to 
dynamic aeroelastic effects. This experiment was conducted in a fixed-base simulation, 
treated the longitudinal axis only and used a pursuit display in which the effects of elastic 
deformation were clearly indicated. 

In the experiment described herein, both the pitch and roll vehicle axes are being 
controlled and the effects of motion cues and display dynamics are specifically addressed. 
The experiment addresses these effects in a generic way for a representative flexible aircraft. 
It is intended that the results help further identify important factors in determining overall 
aeroelastic effects. 

The experimental results include complete response data, as well as subjective data, to 
allow a variety of analyses. The data was recorded in several forms to facilitate subsequent 
analyses. These include vehicle responses recorded on magnetic tape, and pilot comments and 
subjective ratings with complete cross-referencing of all the data. Simulator frequency 
response data was also collected to more completely document the simulation facility and 
provide a database that may be used later to address simulation fidelity issues. 

The simulation facility that was used is the NASA Langley Visual/Motion Simulator 
(VMS). This simulator has programmable displays and control manipulators along with 
motion capability. The elastic aircraft model that was used was developed by modifying an 
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CHAPTER 2 
THE EXPERIMENT 

The documentation of the experiment involves descriptions of - 1) the simulation facility; 
2) the vehicle math model; and 3) the tasks, displays and configurations. 

Facility 

The NASA Langley Visual/Motion Simulator (VMS) Facility was used for the 
experiment and is described in the functional block diagram shown in Figure 1. Each of the 
blocks represent dynamic elements in the simulation. This structure is useful, not only for 
visualizing the various functions of the simulator, but also for understanding the dynamics of 
the simulator itself. 

The VMS is a general purpose simulation facility consisting of a two-man fully 
instrumented cockpit mounted on a six-degree-of-ffeedom synergistic motion baseJ 5,6,7 ' The 
computational functions of the simulator are accomplished by means of a CYBER 175 
computer The visual display and conventional cockpit instrumentation are capable of 
accepting inputs from several sources. 

The control manipulators that were used include a center stick for pitch and roll control, 
rudder pedals for yaw control, and a throttle. The instrumentation included the standard 
aircraft instruments; airspeed indicator, turn indicator, artificial horizon, vertical speed 
indicator, etc. The command display was implemented with a head-up type of display using a 
CRT. 

The control stick and pedals effect command inputs through a programmable hydraulic 
control loading system. This system provides control of the force/feel feedback to the pilot 
and was used to simulate realistic stick and pedal forces and travel distances. The control 
loader was adjusted to provide good force/feel feedback for the baseline vehicle and was not 
varied among the configurations. The values that were used are presented in Table 1. 

The display of the VMS is driven by the Graphics Display System which consists of two 
Adage Graphics Terminals (AGT’s). Each AGT has three display generators which are 
interfaced with television cameras to provide cockpit displays. Each AGT is also digitally 
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existing NASA Langley simulation model of a rigid vehicle and was implemented to allow 
on-line variation of several configuration parameters. Among these were, for example, 
invacuo-vibration frequencies allowing the experimenter to quickly alter the effective stiffness 
of the vehicle. The tasks and configurations were designed to provide a challenging 
environment for the pilots. 

This report is organized to systematically present the description, documentation and 
results of the simulation experiment described above. First, the simulation facility, vehicle 
model and experimental set-up are discussed. Next, the basic results are presented and 
discussed. In addition, some preliminary analysis of the data is presented along with some 
discussion and interpretation of the results. The report then concludes with a general summary 
and suggestions for future applications and analyses of the data. 




Figure 1 - Simulator Functional Block Diagram 
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Table 1 - Force/Feel System Parameters 




Breakout 


Friction 

Force 




Travel 

Force 

Preload 

Force 

Gradient 

Trim Limits 


(± in.) 

(lbs.) 

(lbs.) 

(lbs.) 

Obs/in) 

Position 

Rate 

Pitch 

6 aft 





2.7" fwd stick 




3.3 

1.8 

1.5 

10 


0.5 — 
sec 

Stick 

4 fwd 





5.3" aft stick 


Roll 

±5 

3.0 

1.4 

1.6 

6.5 

±20° of tail 

2.0^1. 

sec 

Stick 






deflection 


Pedal 

±3.25 

14.0 

6.5 

mm 

40.0 

±1.3" pedal 

0.3 — 
sec 
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interfaced with the CYBER 175 to provide the capability to drive the displays with simulation 
variables during real-time simulation. 


Vehicle Model 

The vehicle that was modeled (shown in Figure 2) is similar to an early version of the 
Rockwell B-l and is considered to be typical of advanced transport- type aircraft. This vehicle 
was chosen also because of the availability of structural data and of a nonlinear rigid-body 
simulation model that was used in a previous NASA study. 


Airframe 

The rigid vehicle model consists of the rigid-body equations of motion, tabulated 
nonlinear aerodynamic coefficients, control system and engine models, auxiliary and 
kinematic equations, and force/feel feedback (i.e. control loader) parameter values. The rigid 
model was then appended with structural dynamics which are coupled to the rigid-body 
degrees-of-freedom via the aerodynamic forces. 

The equations of motion for the rigid model are written in a body-fixed reference frame. 
They consist of the traditional six rigid-body degrees-of-freedom, three translational and three 
rotational. The structural dynamics that were added to the model consist of two structural 
degrees-of-freedom which correspond to a symmetric and antisymmetric free vibration mode 
of the structure. 

The integration of the rigid-body degrees-of-freedom with the aeroelastic degrees-of- 
freedom was accomplished according to the modeling method discussed in Reference 9. The 
equations which result describe the overall motion of the vehicle relative to a mean-axes 
coordinate system, with origin at the instantaneous center of mass. The structural 
displacements are assumed to be small so that variations in the moments of inertia are 
negligible. 




Figure 2 - Geometry of Study Vehicle 
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The equations of motion for the simulation model of the flexible vehicle are presented in 
Table 2. Notice that the form of the equations is very similar to that of a rigid aircraft 10 ^ but 
with the addition of equations associated with structural degrees-of-freedom. This allows the 
rigid vehicle model data to be used with very little alteration. The primary difference between 
the rigid model and the elastic model, beside the additional equations, is in the aerodynamic 
forces and moments (viz. X, Y, Z, L, M, and N). The coupling between the rigid-body 
degrees-of-freedom and the aeroelastic degrees-of-freedom occurs via the aerodynamic forces 
and moments. 

Consider, as an example, a flexible vehicle model for the longitudinal degrees-of- 
freedom with one structural mode in straight, horizontal flight. The equations of motion for 
this model are, 

U = -qW-gsin9 + X (la) 

W = Uq + gcos(|)cos0 + Z (lb) 


q = M 


(lc) 


r) = -2 £cot| - co 2 ri + 


(Id) 


where U and W are the surge and plunge velocities in the body-reference axes, respectively; q 
is the pitch rate; <(> and 9 are the Euler angles of the vehicle; T|, £, and co are the generalized 
coordinate, damping and frequency for the structural mode, respectively; and X, Z, M and 
are the aerodynamic forces, moment and generalized force, respectively. Assuming the 
vehicle has elevator deflection as the only control variable, the aerodynamic force, Z, for 
example, might be expressed as 


pv 0 2 s r 




2M L 


C Zo +C Z« a + C Z s 5 + 



pVpSc 

4M 


+ + c z^ T l 


( 2 ) 


where is the aerodynamic stability derivative associated with the variable (•), p is the air 
density, V 0 is the vehicle velocity, S is the planform area, c is the mean aerodynamic chord, 
and M is the mass of the vehicle. The coupling between the aeroelastic degree-of-ffeedom 
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Table 2 - Vehicle Equations of Motion 

U = r V - q W - g sin 0 + X 
V = pW - r U + g sin <|> cos 0 + Y 
W = Uq-Vp + g cos <j> cos 0 + Z 

P=y- {(Iyy-^qr +Ixz(r + Pq)}+t 

4 = J- ta a -I,x)pr +I«(r 2 -p J »+M 

lyy 

r = ~ {I„ - + Wp- q r » + N 

A zz 

T|y = — 2£yC0y'n y ~ (O^T|y 4" QfJ y 

r\ z = -2C 2 to z ri z - co z 2 n z + Q n , 
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and rigid- body vertical translation occurs through the coefficients C^, and C^. Similar 

coupling terms arise in the (generalized) forces for each degree-of-freedom, rigid and elastic, 
in the model. Some of these terms, furthermore, depend on the vibration mode shapes (see 
Reference 7). 

The original, rigid vehicle model has several control effectors. These include horizontal 
tail deflection, 5 t , differential horizontal tail deflection, 8 DH , wing spoiler deflection, 5 S , upper 
split-rudder deflection, 5 RU , lower split-rudder deflection, 8 RL , and thrust, T (via power lever 
angle, PLA). Additionally, estimates for the control power available from the forward control 
vane, 8 CV , were included in this study. 

The expressions for the aerodynamic forces and moments and generalized forces are 
shown in Table 3. The numerical values for the aerodynamic force coefficients are presented 
in Appendix 1. This tabular data was used with a “table look-up scheme” in the simulation to 
determine the aerodynamic forces and moments. 

The geometry, mass and inertia, and flight condition for the simulation model are shown 
in Table 4. The shapes of the two structural modes included in the model are shown in Figure 
3. More complete modal data for the baseline vehicle is presented in Appendix 2. This 
includes finite element control point locations and control point displacements for the two 
modeled structural modes. 

Additional information needed to drive the displays, instrumentation and motion base of 
the simulator are given in Table 5. 


Augmentation 

The simulation model includes actuator dynamics and stability and control augmentation 
system (SCAS) dynamics as well. These were necessary to assure good handling qualities for 
the baseline vehicle. Block diagrams of the SCAS systems for the pitch, roll, and yaw axes 
are shown in Figures 4 through 6. The actuator/servo elements that are modeled consist of 
lags, limiters, and gearing effects. The control laws include both constant and variable gains, 
scheduled with altitude and Mach number. The constant gain values are tabulated in Table 6 
and the scheduled gains are presented in Appendix 3. The limiters and gearing effects are also 
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Table 3 - Aerodynamic Forces and Moments 


pv 0 2 s 


pV 0 Sc 


[ C X+ C XsA+Cx 8. 1 5 s l+ Cx8 cv S cv ] + [ C XqQ] + 


Y = -^^[Cyp+Cys^SrkKTys^Srl+Cys Ss+Cys^bDirHCY^^y] 


pV 0 Sb 


[C Y ^y] 


PV 0 2 S 


pV e Sc 

[Cz+C z5 8 t +C z5i 1 8 S I +C Z 5 c< 8 cv +C Zlli r| J + — [Czqq+ c zn, T lJ 


pV,?Sb 

L = — — [C / p(3fC /SRU 5 RU +C /8RL 8 RL +C/8 8 s +C/ 5DH 8 DH +C /Tly , n y ] 


pV 0 Sb 2 


■[C /p p+C /r r+C /T1 Tiy] 


pV 2 Sc 


2Iyy ^ ^ 8 g I +CmS cv ^cv^’ ^-■ m T) 1 ^lz] 


pV„Sc 2 


+ [Cm,q+c m „«^.nj + ^ t 


pV 0 2 Sb 


[Qlf P + Qi5|Ru^RU^^5iu.^RL^n8 > ^s^n5 D H^Dir*^Trny 1 1y] 


pV Q Sb 2 . 

+ ^^^[CnpP^+CHp^C^y] 


pV 0 Sb 

+ Q|y8 c »^CV + QlyTly r ly] + 2Af y ^>pP 4 ^>^^lxT|y T *y^ 

Qn» = 2M 

pV-Sc 

+ 4 M z ^ C ’i^ q+C ^ a+C n.n, 1 lJ 
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Geometry 


Weight 


Modal 

Generalized Mass 


Flight 

Condition 







1.34 ft. (mean chord) 
136.7 ft. (wing span) 
1946 ft 2 (planform area) 
65° (sweep angle) 


288,000 lbs (gross weight) 


950.000 slug-ft 2 

6.400.000 slug-ft 2 

7.100.000 slug-ft 2 
-52,700 slug-ft 2 


28,991 slug-ft 2 (antisymmetric model) 
183.6 slug-ft 2 (symmetric mode) 


0.002055 (air density) 
ft 

657.5 (cruise velocity) 

sec 


















Figure 3 - Shapes of Modeled Structural Modes 





Tabic 5 - Additional Vehicle Equations 


Euler Angles and Auxiliary Equations 
0 = q cos $ - r sin 

<|>=p+ q tan 0 sin $ + r tan 0 cos <)> 

Y = r cos $ sec 0 + q sin $ sec 0 

V 0 = [U 2 +V 2 +W 2 ]* (velocity) 
a = Tan -1 (-^-) (angle of attack) 

P = Tan -1 (-^-) (side slip) 

*0 

h = U sin 0 - V cos 0 sin <{> - W cos 0 cos <f> (rate of climb) 


Acclerations at Center of Mass (g’s) 


a. = — (U - V r + Wq + g sin 0) 


8v = — (V - Wp + Ur - g cos 0 sin $) 

g 


a^ = y (g cos 0 cos <> - W + Uq + Vr) = -a^ 




Accelerations tt Cockpit (g’s) 
- (q + Pr) - (r 2 + q 2 ) — 


•y, = a y, + ^ r + P < 1) ~- + (P - r 9 + ”?■ 

6 o 




<D>p = 0.45 (ft) 
<!>£, = 0.39 (ft) 


Numerical Values at Cockpit 

icp = 0.0 (ft) 


lcp*63.4 (ft) 


<!>♦, = -0.0027 (rad) 




Primary 
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Figure 4 - Pitch SCAS Block Diagram 
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Figure 5 - Roll SCAS Block Diagram 
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Figure 6 - Yaw SCAS Block Diagram 
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Table 6 - Constant Control Gains 


Pitch 

SCAS 

v c c DEG 
Kt 55 IN 
Xj = 0.02 sec 

Kq = 1.6 or 2.0 sec (configuration dependent) 

Kq = 4.58 DEG 
g 

Roll 

SCAS 

Ko„ = 9.0 DEG 

K 2 = 0.02 sec 

K — ft 57 

^ * 7 deg/sec 

Yaw 

K, v = 200.0 

SCAS 

, g 

K, = -7.7 
^ in. 

Phugoid 

Augmentation 

K pa = 250 "7“ 

VA ft 
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presented in Appendix 3. It is important to note that for the most part, the control system 
gains were not variables of the experiment. The control gains were altered on very few 
configurations and are so noted. 

An additional controller was added to the flexible aircraft model to increase the damping 
of the phugoid mode for the baseline configuration. This was done to improve the handling 
qualities of the baseline configuration and reduce the need for the pilot to monitor speed. The 
phugoid augmentation controller consists of a velocity error signal and gain, as shown in the 
block diagram in Figure 7. The value of the gain, K PA , is shown in Table 6 with the other 
controller gains. Note that, in the simulation data this phugoid augmentation is referred to as a 
“speed hold” controller. 

To illustrate the effects of varying the structural stiffness (i.e. changing the invacuo 
vibration frequencies) on the vehicle dynamics, pitch-rate-to-stick-deflection transfer functions 
are shown in Table 7, both for total pitch rate measured at the cockpit location, q T , and for the 
pitch rate of the ideal body-reference (mean) axes, qj^. These transfer functions were obtained 
from the longitudinal degrees-of-freedom of the baseline nonlinear simulation model, 
linearized in level flight at the flight conditions summarized in Table 4. The transfer functions 
include the effect of the baseline pitch SCAS, phugoid augmentation and a first-order 
approximation for the actuator lag, with a time constant of 0.1 seconds. Engine dynamic 
effects, however, are not included. 

Two of the modes roughly correspond to the ' standard dynamic modes of aircraft 
dynamics. These include a phugoid-like mode and a short-period-like mode. In addition, an 
aeroelastic-structural mode and an (augmented) actuator mode are present. The poles 
associated with these modes are listed in Table 7 in the following order; phugoid mode, short 
period mode, aeroelastic mode and augmented actuator mode. These dynamic modes are 
readily identified for the stiffer configurations, that is those with a symmetric mode frequency 
of 1.5 or 2.0 Hertz. However, as the symmetric mode frequency is reduced the dynamics 
change from these “standard” modes. In particular, the “short period” mode and 
“structural” mode become highly coupled. Note also that the phugoid poles begin as a well 
damped complex conjugate pair but migrate to become two real poles as the vibration 
frequency is reduced. In fact, one phugoid pole becomes slightly unstable in the most flexible 
configuration. 
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Table 7 



Transfer Function Data for Simulation Model 



Invacuo Symmetric Mode Frequency (Hz) 


2.0 

1.5 

1.0 

0.8 

T.F. Gain 

4.764 

4.764 

4.782 

4.782 


0.0 

0.0 

0.0 

0.00 

q x 

*p 

-0.0647 

-0.0666 

-0.0511 

-0.0652 

Zeros* 

-0.440 

-0.491 




-5.1205 

-3.569 


(1.211,0.744) 


4.499 

2.985 

1.347 

0.578 

T.F. Gain 

-0.858 

-0.858 

-0.839 

-0.839 


0.0 

0.0 

0.0 

0.0 

< 1 m 

-0.0624 

-0.0620 

-0.0361 

-0.0365 

8 p 

Zeros* 

-0.402 

-0.408 

-0.422 

-0.509 


(11.86,0.008) 

(8.45,0.003) 

(4.77,-0.013) 

(2.77,-0.047) 


(0.0397,0.675) 

(0.037,0.748) 

- 0.0010 

-0.0229 

-0.0495 

0.0363 

Poles* 

(6.054,0.820) 

(5.615,0.825) 

(4.03,0.795) 

(2.14,0.453) 


(12.57,0.0544) 

(9.46,0.081) 

(6.523,0.281) 

(6.359,0.460) 


-1.733 

-1.857 

-2.322 

-4.123 


♦complex conjugate poles and zeros written: (to (rad/sec), £) 





































22 


The transfer functions for total pitch rate, q T , are all nonminimum phase due to the mode 
shape of the vibration mode, that is, structural deflections cause initial pitch displacements to 
be opposite in direction to the final pitch attitude. The location of the nonminimum phase zero 
varies considerably as a function of symmetric mode frequency. In addition, the migration of 
the other zeros due to changes in the structural mode frequency also affect the vehicle 
dynamics, especially the two negative real zeros which merge and break off into two well 
damped oscillatory zeros when the symmetric mode frequency is reduced to 0.8 Hertz. 

The set of transfer functions for mean-axis pitch' rate, q M , are fundamentally different. 
They are minimum phase and have a pair of lightly damped oscillatory zeros that for the 
stiffer configurations are located close to the complex poles associated with the aeroelastic 
mode. As the frequency of the symmetric mode is reduced, the oscillatory zeros migrate until, 
for the more flexible configuration, they are closer to the two low-frequency complex poles. 
This clearly reflects fundamentally different dynamics, compared to the stiffer configurations. 

An additional exploratory control system was briefly implemented in the simulation as a 
demonstration. The structure of the simulation model, both in the structure of the equations of 
motion and in the structure of the SCAS implementation allows the model to be readily altered 
by appending additional system dynamics. This demonstration controller was synthesized 
using eigenspace assignment methods^ 1 ^ and implemented in place of the pitch SCAS 
controller. Although the performance of such a control law was investigated briefly, handling 
qualities assessments were not obtained. 

The implementation of the eigenspace assignment controller was accomplished as shown 
in Figure 8. The controller requires feedback of the velocities, U and W, pitch attitude and 
pitch attitude rate, 0 and q, plunge acceleration and plunge acceleration rate, n z and ri z , tail 
deflection, 8 t , and control vane deflection, 8 CV . The command signals for this controller are 
determined from Equations (3) and (4). The control gains for each of the parameters, given in 
these equations, are presented in Table 8. 


Non-Linear 
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Figure 8 - Eigenspace Assignment Controller Block Diagram 
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Table 8 - Eigenspace Assignment Controller Gains 


Tail Command Gains 

Vane Command Gains 

K‘ 

- 0.51172 

K u cv 

-1.2567 

K 

0.07557 

XT CV 

1V W 

0.1498 

K 

1.5044 

K^ v 

-61.693 


-85.842 

K« cv 

-192.59 

K 

0.4197 

K 

4.1918 


0.1336 

K iT 

0.3579 


-0.1575 


0.4970 


-0.02579 


-0.09404 


• • • **Q/1 

(note: all 8’s in degrees; u, w in ft/sec; 0,r| in ; 0, r\ in rad) 

sec 
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Sta» = S.™, - K'Ju - U TRIM 1 + K‘ w [ W - W TRIM ' 

k. J L J 


+ K l n 


Q — Qtrim + K e 


+ K\ 


0 “ ®TRIM 

I 

"z “ ^zTRIm] 


+ K l ; 


6, 


^CV ^cv-nuM 


i] +K l , 

"z"^™,] + K i. 

VO +Kt 6„ 

. J 

6cv com = 5 CVtom - K- u [u - U TRIM ] + K~ w [w - W TRIM ' 
+ K cv q |^q - qxRiMj + K cv e ^0 - 0 trim 

+ K< \ [ n z“ n *nuii] + [ n z - n zTRIM 

^[5 t -5 tTOM ] + K CV 5 


+ K CV 


^cv ^cvtwm 


(3) 


(4) 


Engine 

The simulation model also includes engine dynamics. This model was present in the 
original, rigid model and was retained in the flexible model. A block diagram is shown in 
Figure 9. The functional elements (e.g. f(PLA,h,M) and f(h,M) ) which appear in the model 
are expressed in the tabular forms presented in Tables 9 through 11. It should be noted that 

I 

the effect of the engine dynamics on the simulation results are most likely very small, and the 
throttle was utilized little by the pilots during the experiment. 


Tasks and Displays 


The multi-axis piloting tasks used in the simulation experiment were developed to 
require the pilot to fully utilize his entire dynamic range, to provide an indication of the pilot’s 
ability to excite the aeroelastic modes. The tasks designed were generic, but intended to be 
representative of activities experienced in actual flight. The intent is to allow the results of the 

simulation to be better interpreted in terms of actual vehicle requirements and to allow the 
pilot to relate the results to his previous experiences. 


36089 
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Table 9 - Engine Data : FNE 


^1 


.3 

.4 

.5 

.6 

.7 

.8 

.9 

0 

26000 

28400 

29300 

30500 

31700 

32150 

32200 

31900 


13800 

13800 

13800 

15300 

16800 

18500 

20300 

22200 

36089 

7900 

7900 

7900 

7900 

7900 

8900 

10000 

11100 


Table 10 - Engine Data : %FNE / sec 


%FNE 

THCINC 

%FNE 

THCDEC 

-8 

2 

-8.2 

-1.3 

-3 

2.63 

-5.5 

-1.8 

1.5 

7.5 

.5 

-12 

11.5 

10.53 

10 

-19 

14 

2.38 

16.5 

-6.5 

19.5 

5.5 

36.5 

-20 

31 

11.5 

58 

-43 

100 

19.71 

90 

-60 

200 

.01 

95 

-33.33 

300 

.01 

100 

-14.29 













28 














































29 


The two generic tasks developed to fulfill these objectives are 

(1) a pitch and roll attitude pursuit tracking task, and 

(2) a flight-path elevation and heading tracking task. 

These tasks are representative of the critical loop closures in the pilot/vehicle system for many 
aircraft operations.^ The pitch/roll task was chosen to evaluate characteristics of the 
“inner-loop” dynamics and the flight-path/heading task was chosen to evaluate characteristics 
of the overall, “outer-loop” dynamics of the vehicle. The complement of piloting tasks that 
were utilized during the simulation experiment primarily included these two dual-axis tasks, 
along with some single-axis variants of these tasks. 

The visual displays that were used in the tasks were designed to provide the pilot with the 
necessary information in a way that was complete yet made use of simple symbology. In 
addition, the displays had to be acceptable to the pilot and provide the information in a 
consistent, “realistic” manner. 

In order to accomplish these objectives, the displays were chosen to be of the inside-out 
and “fly-to” format. That is, the own-ship symbol is fixed on the display, the horizon is 
driven by the aircraft pitch attitude and bank angle, and the target is driven by the tracking 
errors. The pilot attempts to fly to the target. 

The symbology that was used in the displays is similar to that depicted in Figure 10. 
This symbology was implemented via the AGT and displayed to the pilot by means of a CRT 
in the simulator cockpit. The display parameters that determine the display motion scaling, 
target range, field of view, etc. were chosen to provide sufficient resolution of attitude and 
error information for the task. The display scaling used is presented in Table 12. 

The pitch/roll tasks are pursuit tasks. That is, both pitch and/or roll attitude errors and the 
pitch and/or roll attitudes of the aircraft is displayed to the pilot. The horizon is driven by the 
displayed pitch and roll attitude responses, 0 D and <{> D , but the vehicle responses used were 
varied. A discussion of these variations will be presented later. The target aircraft relative roll 
attitude and position are described by the pitch and/or roll tracking errors, ee and e^, where 
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Figure 10 - Tracking Task Display Symbology 
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Table 12 • 

■ Display Scaling Parameters 

Parameter 

Scaling 

®D* Eq. 

Displacement of horizon and 
target yields one-to-one 
angular deflection at pilot’s eye 


Angular rotation of horizon 
on display is one-to-one 

«♦ 

Angular rotation of target 
on display is two-to-one 


Table 13 - Command Parameters 



Tracking Task 

Parameter 

Pitch/Roll 

Flight-Path/Heading 

t\j (mean) 

Orad 

Orad 

ri 2 (mean) 

Orad 

Orad 

o 2 t (variance) 

0.4/32 rad 2 

0.3/32 rad 2 

(variance) 

0.77/32 rad 2 

1.5/32 rad 2 


0.5 rad/sec 

0.3 rad/sec 

0)$ 

0.5 rad/sec 

0.3 rad/sec 


0.5 sec 

0.5 sec 

a e c 

2.7 deg 

1.8 deg 

°«c 

3.8 deg 

4.1 deg 
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£q = 0 C - 0 D 

(5a) 

e <j> = 0C - 0D 

(5b) 


and where 0 C and <J) C are the command signals. 

The command signals that were used in the task were chosen to provide the pilot with a 
challenging task, with the potential to result in excitation of the aeroelastic modes. They were 
driven by two independent random signals, governed by the following relations. 

0 C + 1.4to e 0 c + (o 0 2 e c = o) e 2 r| 1 (6a) 

ic + 1 -4 co^ <j> c + V 0c = O)^ 2 r| 2 (6b) 

The statistics for the random signals, % and TI 2 , were chosen to obtain command signals with 
the desired statistics. The appropriate statistical parameters and filter frequencies are 
presented in Table 13. 

The flight-path and heading tracking tasks also provided the pilot with the aircraft 
attitude situation as well as the flight-path tracking errors. The displayed attitudes (horizon) 
for the pilot’s vehicle in these tasks are once again 0 D and 0 D . But, the target aircraft relative 
positions are described by flight-path and heading errors, £y and e^, where 



ty = [ Yc “ y] cos0 D + [ v c - \j/j sin<j) D 

(7a) 

and where, 

£V= [ [ Yc-v) — By sin \j/ D J /cos0 D 

(7b) 


y=0 D -a 

(8a) 


I 

> 

II 

l> 

(8b) 


where \j/ 0 is a reference heading. The command signals for flight-path and heading, y c and 
Y c , are related to the pitch and roll commands through the kinematic relations shown below. 
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Y=-±y+±e c 

x c x c 
Yc = / cos <{> c 


(9a) 

(9b) 

(9c) 


¥c = Y' sin<])c 

The appropriate values for x c are also given in Table 13. 

The displayed attitudes of the aircraft, 0 D and <J) D , were variable as part of the 
experiment, leading to two different types of display dynamics. The first represents an “out- 
of- the- window” display, that is the displayed attitudes correspond to those of the aircraft at 
the cockpit location, including the effects of structural displacements. In this case, 

0 D = 0 + KeTi z (10a) 

<t>D = ( t ) + K <|) T ly (10b) 

where 0 and <J) are the rigid-body (mean-axis) Euler angles of the aircraft, K 0 and are the 
slopes of the two modeled structural modes evaluated at the cockpit location, (i.e. Kq = -0.03 
rad, and = -0.0027 rad), and r) y and tj z are the generalized coordinates of the two modeled 
structural modes. 

The second variation can be thought of as a filtered display. In this case, 

0 D = 0 (11a) 


(t) D -<t) 


(Hb) 


or the dynamic structural effects are filtered out. In subsequent discussions regarding the two 
display variations described above, the former (Equation 10) will be referred to as the 
“flexible display” and the later (Equation 11) will be referred to the “rigid display.” 
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Configurations 

The experimental configurations that were used in the simulation study were developed 
to address variations in structural stiffness, simulator motion effects, display dynamics effects, 
and consequences of multi-axis versus single-axis tasks. As mentioned previously, the 
variation in structural stiffness was accomplished by changing the invacuo free vibration 
frequencies of the two modeled structural modes. The baseline vehicle had invacuo vibration 
frequencies of 2.0 Hertz for each mode. These values were varied between 0.8 and 2.0 Hertz 
for the symmetric mode and between 1.0 and 2.0 Hertz for the antisymmetric mode. The 
simulator motion effects were examined by conducting some simulation runs with the same 
model dynamics, both with the motion enabled and disabled. Similarly, the effects of display 
dynamics were addressed by conducting runs with the same model dynamics, but with both 
display variations. The multi-axis/single-axis task effects were examined by simulating the 
same dynamics in both single-axis and dual-axis tasks. 

The configuration variations described above can be categorized by the choice of several 
parameters. Some of these parameters are associated with the experimental “set-up” (e.g. 
simulator motion, pilot task, etc.) and others are related to the “dynamic properties” of the 
vehicle (e.g. invacuo vibration frequencies). Case numbers were assigned to each 
configuration to provide a systematic way to differentiate between various combinations of 
parameters. 

The set-up parameters are reflected in the first four entries in the case number. These 
“identifiers” are defined in Table 14. An example is also included to demonstrate the 
numbering system. Notice that the last entry in the case number is called a “dynamic 
parameters identifier.” This identifier is used to indicate the combination of dynamic 
parameters chosen for each configuration. Figure 1 1 is an excerpt of the list of configurations 
tested, the entire list is presented in Appendix 4. The combination of dynamic parameters for 
each configuration is defined in this list. 

Notice that in addition to the frequencies of the invacuo vibration modes, the dynamic 
information includes structural mode status (e.g. ETAY OFF implies the dynamics were that 
of the rigid aircraft in the lateral axis), SCAS status, and speed hold (i.e. phugoid 
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Table 14 - Case Number Definition 


Dynamic Parameters Identifier (see Figure 1 1 and Appendix 4) 
Display Identifier f 1 - flexible display 
1 2 - rigid display 
Motion Identifier f 1 - motion on 
L 2 - motion off 
' A - longitudinal 

Task Identifier j B - lateral-directional axis only 
Getter) C - two axes (both long. & lat.-dir.) 

Task Identifier f 1 - pitch/roll tracking 
(number) i 2 - heading/flight path tracking 

Example: Case No. 1C- 12-5 

task properties -pitch/roll tracking task 

-long, (pitch) axis and lat. (roll) axis enabled 
-motion on 

-standard display in use (flexible display) 

dynamic properties -C 0 y = 2.0 Hz (antisymmetric mode freq.) 

-co z = 1.0 Hz (symmetric mode freq.) 
-SCAS-on (Kq = 2.0) 

-Speed Hold-on 
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CASE NO. 

ETAY 

ETAZ 

OMEGY 

OMEGZ 

SCAS 

KQ 

SPDH 

OLD CASE 

1 A— 11—1 

OFF 

OFF 

— 

— 

ON 

2.0 

OFF 

4 

1A-11-2 

OFF 

ON 

2.0 

2.0 

ON 

2.0 

OFF 

10 

1A-11-3 

OFF 

ON 

2.0 

1.75 

ON 

2.0 

OFF 

16 

1A-11-4 

OFF 

ON 

2.0 

1.5 

ON 

2.0 

ON 

22 

1A-21-1 

OFF 

OFF 

— 

— 

ON 

1.6 

OFF 

1* 

1A-21-2 

OFF 

OFF 

- 

- 

ON 

2.0 

OFF 

1 

1 A— 2 1 —3 

OFF 

ON 

2.0 

2.0 

ON 

2.0 

OFF 

7 

1A-21-4 

OFF 

ON 

2.0 

1.5 

ON 

2.0 

ON 

21 

IB-11-1 

OFF 

OFF 

— 

— 

ON 

2.0 

OFF 

5 

IB-11-2 

ON 

OFF 

2.0 

2.0 

ON 

2.0 

OFF 

11 

IB-11-3 

ON 

OFF 

1.5 

2.0 

ON 

2.0 

ON 

32 

IB— 21— 1 

OFF 

OFF 

_ 

— 

ON 

1.6 

OFF 

2* 

IB-21-2 

OFF 

OFF 

- 

- 

ON 

2.0 

OFF 

2 

IB-21-3 

ON 

OFF 

2.0 

2.0 

ON 

2.0 

OFF 

8 

IB-21-4 

ON 

OFF 

1.5 

2.0 

ON 

2.0 

ON 

31 


Figure 1 1 - Configuration Listing (excerpt) 


DATE 

| RUN j 

| CASE | 

| COMMENT 

4/30 

1 

1A-21-1 

no taped data, limited strips 


2 

IB-21-1 

* * 


3 

1C-21-1 



4 

1A-21-2 



5 

1A-21-2 

* + 

5/2 

1 

1A-21-1 



2 

IB-21-2 



3 

1C-21-2 



4 

1 A- 11-1 



5 

1 B— 1 1—1 



6 

lC-11-1 


5/3 

1 

lC-11-1 

tau(3) = 0 


2 

1C— 1 1—1 



3 

1A-21-3 



4 

IB-21-3 



5 

1C-21-3 



6 

1A-11-2 



7 

IB-11-2 



8 

1C— 1 1 —2 



Figure 12 - Chronological Listing of Runs (excerpt) 
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augmentation) controller status. Notice also that the gain on pitch-rate feedback in the pitch 
SCAS, KQ, is indicated. This parameter was varied as part of the tuning process which will 
be discussed in the next chapter. The entry for pitch SCAS gain was also used to indicate 
when the eigenspace assignment controller was used. When this is the case, ESC appears 
rather than a numerical value for KQ. The next item is the the entry entitled “SPDH” (for 
speed hold) and is used to indicate the status of the phugoid augmentation controller. Finally, 
the parameter called “old case” is tabulated. This number corresponds to the case number 
originally assigned to the configuration during the study. It is included in this list since it is 
useful in cross-referencing the experimental data, for example in relating the data recorded on 
the strip charts and the pilot comments made during the simulation sessions. 

To further facilitate cross-referencing of the experimental data, a chronological listing of 
the simulation runs and case numbers has been developed as well. This list presents the case 
number corresponding to each simulation run, it also includes comments regarding the 
availability of data associated with particular runs. An excerpt of this list is shown in Figure 
12 and the entire chronological listing is presented in Appendix 4. 

The baseline vehicle configuration corresponds approximately to the dynamics of an 
early B-l aircraft. In terms of the parameters described above, this chosen baseline 
configuration includes the dynamics of the two structural modes with frequencies of 2.0 Hertz, 
simulator motion disabled, and the flexible display. A more complete specification of the 
baseline configuration is shown in Table 15. 
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Table 15 - Baseline Configuration Parameters 


Flight 

Condition 

A = 65° (swing sweep) 
h = 5000 ft (altitude) 

M = 0.6 (Mach number) 

Vehicle 

Dynamics 

Flexible modes - enabled 
CDy » 2.0 Hz (vibration 

0) z = 2.0 Hz frequencies) 


SCAS - enabled 

Kq as 2.0 sec (pitch rate gain) 
Phugoid Augmentation - enabled 

Motion 
and Display 

Motion - disabled 

Display - flexible ("out-of-window") 

Task 

and 

Commands 

Task - pitch/roll tracking 

Ge = 0.4 (command 
= 0.77 variances) 

x c = 0.5 (sec) 











39 


CHAPTER 3 

RESULTS AND ANALYSIS 


Conduct of the Experiment 

The vehicle configurations were evaluated in the order specified in the chronological 
listing shown in Appendix 4. Three pilots participated in evaluating the experimental 
configurations, although pilot B was the primary test pilot. The background and experience of 
each pilot is summarized in Table 16. Their backgrounds reflect a wide range of piloting 
experience and this should be considered when the simulation results are reviewed. 

Each pilot was briefed on the experimental structure and task objectives. They were 
required to keep the center of the target (i.e. the intersection of the wing-line and the tail-line) 
inside the circle (or dot) at the center of the own-ship symbol. See Figure 10 for an illustration 
of the display symbology. The pilots were encouraged to use whatever technique they 
considered useful in minimizing the displayed error and were allowed to practice with a 
configuration until they felt they were ready to attempt an evaluation run. When an evaluation 
run was completed, the pilot was asked to comment on any aspect of the run he felt was 
noteworthy. For example, he might note his preferences relative to the vehicle responses, 
performance, and the techniques he employed to perform the task. Then he was requested to 
rate the configuration using the Cooper-Harper rating scale. 

Early in the experiment, the simulation was tuned to improve the “feel” of the 
simulation of the baseline configuration. Several model parameters were varied to obtain a 
baseline configuration which had acceptable and “realistic” responses, when evaluated by the 
pilots. This tuning process resulted in the baseline parameters noted in Table 15. During the 
tuning process small variations in some simulation model parameters occurred. When these 
parameters differ from those of the baseline configuration it is noted in the data records. 

Three types of data were collected during each simulation run: (1) pilot comments and 
subjective ratings (Cooper-Harper), (2) strip chart recordings, and (3) vehicle response time- 
series data (i.e. recorded on magnetic tape). The pilot comments and ratings for each 
simulation run were documented using a summary sheet. Figure 13 depicts a typical summary 
sheet and the entire collection of summary sheets is presented in Appendix 5. The sheets are 
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Pilot 



Background 

Student general aviation pilot. Limited ~ 
experience in piloting aircraft, no experience 
with large aircraft. Familiar with experi- 
mental set-up, model characteristics and 
project goals and expectations. No previous 
experience in rating aircraft handling 
qual ities. 


Experienced test pilot. Extensive background 
in flying large aircraft. Basic (general) 
understanding of experimental set-up and 
objectives. Substantial background in 
rating aircraft handling qualities using 
the Cooper-Harper scale. 


Experienced general aviation pilot. No 
experience with large aircraft or with 
rating handling qualities. Basic under- 
standing of experimental set-up and objectives. 
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SPECS : 





V : ON 

V 

CO 

= 2. 0 

SCAS : ON 

K =2.0 

* 

rt : on 
* 

<0 

ft 

II 

to 

• 

O 

SPEED HOLD 

; ON 

OLD CASE 

NO. : 

13 


SCORES s 




COOPER- 





HARPER 

€ =■■ 

lorta 

. 079 

(7 =1.17 

lont 

RATING 

£ =: 
lot 

• 

10 

O =1.14 

tat 

S 


COMMENTS s 


- PILOT CAN CONTROL THE SITUATION AND RIDE 
IS MUCH SMOOTHER 


VISUAL CUES AID IN CONTROL OF THE 
SITUATION 


QUALITATIVELY DISCONCERTING THAT SCORES 
ARE NOT HIGHER (WORSE) SINCE NOTED MORE 
ERROR 


THE PILOT PERCIEVED WORSE PERFORMANCE 
AND HIGH WORK LOAD 


- NOTE : NO DIGITAL DATA 


Figure 13 - Sample Summary Sheet 
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ordered by case number to make comparisons between and among various cases easier. Each 
sheet is cross-referenced with the strip charts and time response data by means of the run 
number. The run number indicates both the day and order in which the run was made. For 
example, RUN 5-10-3 is the third run completed on May tenth. Other information found on 
the summary sheets are the dynamic parameters described in the previous chapter and are 
included to facilitate later analysis. The pilot identifier is also included to indicate which of 
the evaluation pilots participated in each run. These summary sheets are also used to indicate 
variations from the baseline configuration due to the tuning process described earlier. 

During the simulation runs, the strip charts provided a real-time graphical record of 
twenty-four different responses of interest. Table 17 lists the twenty-four parameters that were 
plotted on three recorders. The strip charts are organized chronologically by run number to 
facilitate cross-referencing with the other data. Figure 14 depicts a sample strip chart plot, 
notice that the parameter names, scales and units are noted directly on each response. 

The time histories recorded on magnetic tape consists of the twenty-six responses listed 
in Table 18. The sampling period was 3/32 second and the average run length was about 120 
seconds. A header/title for each data file on the magnetic tapes indicates the run number (i.e. 
date and chronological order) of each simulation run. This header allows the contents of any 
data file to be readily located, identified and accessed, and facilitates cross-referencing with 
the other forms of data. 

The subjective pilot ratings, of course, provide a measure of the handling qualities of the 
various experimental configurations. In addition, several of the recorded responses provide 
other indications of vehicle performance. For instance, the pitch (or flight-path) and roll (or 
heading) tracking errors provide a measure of tracking performance. Stick displacements and 
rates provide measures of pilot activity and workload. Normal and lateral accelerations 
provide an indication of the ride quality and the pilot induced vibration environment. Root 
mean squared (RMS) values for these vehicle responses were calculated for each simulation 
run. 

The pilot ratings and RMS values of the responses discussed above were summarized 
graphically to allow direct interpretation of the experimental results, and these results are 
grouped by configuration class. A configuration “class” is made up of the experimental 
configurations which possess a specified set of characteristics. These include the task 
definition, motion status, display dynamics, and varying degrees of structural flexibility. The 
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Figure 14 - Sample Strip Chart Plot 



Tabic 18 - Tabulated Response Parameters 
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specific characteristics associated with each configuration class are indicated by the 
configuration numbers defined in Table 14 and the configuration listings in Appendix 4. 


Experimental Results 

An illustration of the simulation results is shown in Figure 15. This plot represents the 
Cooper- Harper pilot ratings that were given to each configuration in the class 1C- 11 with the 
level of flexibility characterized by variations in the symmetric mode frequency (i.e. the 
antisymmetric mode frequency is fixed at 2.0 Hertz). Configuration class 1C- 11 consists of 
experimental configurations in which the dual-axis pitch/roll pursuit tracking task was 
performed with the simulator motion enabled and the flexible display dynamics utilized. Each 
plot indicates which configuration class is being considered and which of the two structural 
modes is used to characterize variations in structural flexibility. Plotted along the axis labeled 
“Level of Flexibility” is simply the invacuo free vibration frequency of the mode which is 
varied. When multiple data points are available for a given level of flexibility, all points 
associated with that frequency are connected by a line. 

The effects of variations in structural flexibility can be addressed by evaluating the trends 
or variations which appear on each individual plot. The effects of variations in the other 
experimental variables, (e.g. task definition, motion status and display status), can be 
addressed by comparing and contrasting the plots for the appropriate configuration classes. 
For example, to compare the effects of motion on the dual-axis pitch/roll tracking task one 
might compare the results for configuration classes 1C- 1 1 and 1C-21 (i.e. those utilizing the 
flexible display dynamics) or configuration classes 1C-12 and 1C-22 (i.e. those utilizing the 
rigid display). 

The complete graphical results are presented in Figures 16 through 26. Figures 16 
through 18 depict the plots of the pilot ratings that were collected during the experiment. 
Figure 16 depicts the pilot ratings associated with the pitch/roll pursuit tracking task with the 
symmetric mode frequency varied. Figure 17 depicts the ratings associated with the same 
task, but with the frequency of the antisymmetric mode varied. Figure 18 depicts the ratings 
associated with the flight-path/heading task where the symmetric mode frequency is varied. 
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Figure 15 - Cooper-Harper Rating : 1C-11 
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Figure 17 - Cooper-Harper Rating : 1C-11, 1C-12, 1C-21 and 1C-22 (Antisymmetric Mode Varied) 
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The ratings are based on the Cooper-Harper scale, and each of the three evaluation pilots is 
indicated with a unique symbol. It should be noted that none of the pilots evaluated every 
configuration. Pilot B was able to evaluate most of the configurations, but pilots A and C only 
evaluated a portion of the possible configurations. Note that multiple data points with the 
same pilot ratings are indicated by multiple rating symbols. Note also that there are many runs 
in which no pilot ratings are available, this is due to the fact that the pilot was only requested 
to rate a given configuration when he felt proficient in the task. However, other data was 
collected for non-rated runs. 

Figures 19 through 26 present the plots of the RMS values for the experimental responses 
described above. Figure 19 presents all the available pitch and flight-path tracking errors. 
Similarly, figures 20 through 26 present the available results for the other responses. The 
RMS pitch (or flight-path) and roll (or heading) errors correspond to the total tracking errors 
(in degrees) and are calculated from the responses P £ and R e defined in Table 18. The RMS 
longitudinal and lateral stick displacements correspond to the linear displacements (in inches) 
of the control stick in the appropriate directions and are calculated from the responses 8 long 
and defined in Table 18. Similarly, the RMS longitudinal and lateral stick rates 
correspond to the linear displacement rates (in inches/second) of the control stick and are 
calculated from the responses 8i ong and 5 lat . Finally, the RMS normal and lateral accelerations 
correspond to the accelerations (in g’s) at the cockpit of the vehicle model (i.e. the 
accelerations that result from the numerical integration of the vehicle equations of motion, not 
the accelerations actually experienced by the pilot which are attenuated due to the washout 
logic and motion-base hardware). These RMS values are calculated from the responses a^ 
and av defined in Table 18. 

J cp 

Note that the RMS values for pitch (or flight-path) error, roll (or heading) error, 
longitudinal stick displacement, lateral stick displacement, longitudinal stick rate, and lateral 
stick rate are plotted for every configuration class for which response data was available. The 
normal accelerations, however, are only presented for variations in symmetric mode frequency 
and the lateral accelerations are presented for variations in antisymmetric mode frequency. 
The letters adjacent to each data point correspond to the pilot identifiers defined in Table 16. 
The asterisk that appears next to some, of the letters indicates that the corresponding 
simulation run was given a subjective pilot rating and should be weighted accordingly. 
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Figure 19c - RMS Pitch and Flight-path Tracking Errors (degrees) 
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Symmetric Mode Frequency - varied Antisymmetric Mode Frequency - 2.0 Hz 
Figure 21d - Longitudinal Stick Displacements (inches) 
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Discussion of Results 

Some comparisons of the graphical results for pilot ratings and RMS responses will be 
considered now. In particular, the results which correspond to the configuration classes 
associated with the dual-axis pitch/roll pursuit tracking task in which the symmetric mode is 
changed to represent variations in structural flexibility will be reviewed. The case number 
designations associated with these configuration classes are 1 C- 11, 1C-12, 1C-21, and 1C- 22. 
Class 1C- 11 corresponds to motion enabled with flexible display dynamics, class 1C- 12 
corresponds to motion enabled but with rigid display dynamics, configuration class 1C-21 
corresponds to motion disabled with flexible display dynamics, and class 1C- 22 corresponds to 
motion disabled with rigid display dynamics. 

Figure 16 depicts the pilot rating results for the specified configuration classes. First look 
at the variations in pilot rating for changes in structural flexibility. A clear trend toward 
increased Cooper-Harper pilot rating (degraded handling qualities) with increased flexibility is 
evident for configuration classes 1C- 11 and 1C-21, those classes with the flexible display 
dynamics. However, the configuration classes which utilize the rigid display (1C- 12 and 1C- 
22), do not exhibit this behavior. In fact, the data for classes 1C- 12 and 1C-22 seem to 
indicate a minimum (best) handling qualities rating for an intermediate level of flexibility (co = 
1.4 Hertz). These trends are indicated regardless of whether motion was used or not. So 
motion appears to have little effect on these pilot ratings, while display dynamics significantly 
altered the pilot’s opinion of the handling qualities. 

Now consider Figure 19b which depicts the RMS pitch tracking errors. The variations in 
tracking errors with changes in structural flexibility indicate the same trends as the pilot 
ratings. Also, these trends are reflected in the roll tracking errors shown in Figure 20b. The 
fact that motion seems to have an insignificant effect on tracking performance is believed to 
be attributable (in part) to the attenuation of the motion commands by the motion logic 
(discussed later). 

Figure 21b depicts the RMS longitudinal stick displacements for the same configuration 
classes (1C-11, 1C-12, 1C-21, and 1C- 22). There appears to be a slight increase in stick 
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displacements with decreased stiffness for all classes. Also, there appears to be little 
difference between flexible and rigid display (1C-X1 versus 1C-X2), or between motion and 
no motion (1C-1X versus 1C-2X). 

Figure 22b depicts the RMS longitudinal stick rates for the specified configuration 
classes. The slight trend toward increased stick rate activity with reduced stiffness is noted 
with rigid display dynamics, (i.e. in configuration classes 1C-12 and 1C-22). On average, the 
more flexible configurations (i.e. with symmetric mode frequency below 1.7 Hertz) have RMS 
stick rates for the rigid display configurations which are clearly higher than those for the 
flexible display configurations. The filtered rigid display may allow the pilot to be more 
aggressive and achieve better performance. There is, however, no apparent effect of motion 
status on longitudinal stick rate. 

Figures 23b and 24b depict the RMS lateral stick displacements and displacement rates 
for the same four configuration classes. Except for a slight reduction in rate with increased 
flexibility for class 1C-11, these responses do not indicate any significant trend between the 
amount of lateral stick activity and the degree of structural flexibility, the presence of motion, 
or the type of display dynamics. Note, however, that structural stiffness for these results is 
characterized by variations in the symmetric mode frequency which should affect the 
longitudinal axis more than the lateral axis. 

A limited spectral analysis of the pilot stick inputs was also performed. The frequency 
spectra, or more specifically the power spectral density (PSD) was calculated from the time 
series data using the Weiner- Khinchine Theorem described in reference 13. Two sets of data 
will be considered. The first set of vehicle configurations serve to address the changes in the 
frequency spectra of pilot inputs due to the effects of motion and display dynamics. The 
second set of vehicle configurations serve to address some relationships between the frequency 
spectra of pilot inputs and pilot ratings, pilot comments and task performance and the effects 
of flexibility. 

First consider some consequences of motion and display dynamics on pilot inputs. For 
this purpose four representative vehicle configurations are compared, each with the same 
dynamics and differing only in motion and display status. Configurations 1C-11-5, 1C-12-3, 
1C-21-5, and 1C-22-2 each have the same dynamics; (pitch/roll) tracking task, symmetric 
mode frequency of 1.5 Hertz, antisymmetric mode frequency of 2.0 Hertz, pitch SCAS pitch- 

c 
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rate gain of 2.0, and phugoid augmentation engaged. They differ however in that 1C-1 1-5 has 
motion enabled and utilizes the flexible display, 1C- 12-3 has motion enabled but utilizes the 
rigid display, configuration 1C-21-5 has the motion disabled and utilizes the flexible display, 
and 1C-22-2 has the motion disabled and uses the rigid display dynamics. The simulation runs 
for each configuration that were used to calculate the PSD’s were performed on the same day 
in consecutive evaluations by the same pilot. 

Figure 27 presents the PSD of the longitudinal stick rate for each of the four 
configurations. Note that for configurations 1C-11-5 (motion) and 1C-21-5 (no motion), those 
configurations that utilize the flexible display dynamics, spectral peaks occur near the 
vibration frequency of the symmetric mode (which is approximately the frequency of the 
associated aeroelastic mode). The other two configurations, with the rigid display dynamics, 
exhibit no peaks of this type. The pilot ratings and tracking performance for the first two 
configurations (flexible display) which demonstrate the spectral peaks had noticeably worse 
dynamic properties. 

The second set of experimental configurations were chosen to help further address 
relationships between the frequency spectra of pilot inputs and pilot ratings, pilot comments, 
and task performance. The configurations used for this analysis include lC-11-1, 1C- 11-2, 
and 1C- 11-3. Each configuration corresponds to the pitch/roll pursuit tracking task. 
Configuration 1C- 11-1 corresponds to a rigid vehicle where the pitch SCAS pitch-rate gain, 
Kq, is 2.0 and the phugoid augmentation is disengaged. Configuration 1C- 11 -2 corresponds to 
a flexible vehicle where both structural mode vibration frequencies are 2.0 Hertz, the pitch 
SCAS pitch-rate gain is 2.0 and the phugoid augmentation is disengaged. Configuration 1C- 
11-3 is identical to configuration 1C- 11-2 with the exception that the phugoid augmentation is 
engaged. Because of their similarity, configurations 1C- 11 -2 and 1C- 11-3 will both be 
referred to as the “flexible configuration” in the subsequent discussion. Configuration 1C- 
11-1 will of course be referred to as the “rigid configuration.” 

Three runs for the two configuration types are considered. Figure 28 presents the PSD of 
the longitudinal stick displacement for each of the three runs for the rigid and flexible 
configurations. This figure also includes the pilot identifier, pilot rating, and RMS pitch 
tracking error for each run. Note that even though the data was obtained from runs that were 
performed on different days, and in one instance by a different pilot, the characteristic shapes 
of the PSD curves are consistent. The rigid configuration has PSD’s which exhibit fairly 
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Stick Rate Configuration: 1C— 12 



Symmetric Mode Frequency - 1.5 Hz Antisymmetric Mode Frequency - 2.0 Hz 


Figure 27 - PSD of Longitudinal Stick Rates - Motion and Display Effects 
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Figure 28 - PSD of Longitudinal Stick Displacements - Rigid and Flexible Configurations 
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smooth, gradual reduction in magnitude as frequency increases. The flexible configuration has 
PSD’s which exhibit higher power below the frequency of the symmetric mode, and a marked 
reduction beyond that point. 

Figure 29 presents the PSD of the longitudinal stick rate for the same simulation runs. 
The pilot identifier, pilot rating, and RMS pitch tracking error are also included here. The 
characteristic shape of the PSD curves are consistent for these responses as well. The rigid 
configuration has power spectra which are fairly constant in magnitude out to a frequency of 
about one Hertz and then drop gradually as frequency increases. The elastic configuration 
exhibits similar behavior but near the symmetric mode frequency (2.0 Hertz) an obvious peak 
occurs which is then followed by a sudden drop in magnitude. 

Now consider the pilot comments and ratings for these configurations in the context of 
the previous discussion. Table 19 presents the pilot comments, subjective ratings and RMS 
pitch tracking errors associated with each of the simulation runs described above. These 
comments indicate that the rigid configuration presented no difficulty in completing the 
required task, consistent with the pilot ratings and tracking errors. However, the flexible 
configuration required the pilot to carefully alter his technique to perform the task 
successfully. Comments such as, 

“let oscillations die out, but could be aggressive,” 

“flying aggressively but smoothly ...” 

“applied learned technique ...” 
and 

“... considerable compensation required” 

indicate that he had to provide fairly complicated compensation to perform the task. The 
comments also tend to imply that the pilot was trying to intentionally avoid exciting the 
structural mode, that is, act like a notch filter. The frequency spectra seem to support this by 
the characteristic shapes of the curves. 

The generally higher pilot ratings that the flexible configuration received (compared to 
those of the rigid configuration) imply that the flexible configuration was significantly harder 
to fly in the specified task. In addition, the variations in pilot rating among the three flexible 
configuration runs (i.e. from Cooper-Harper 4 to 8) indicate that the vehicle is sensitive to the 
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Tabic 19 - Pilot Comments, Ratings and Tracking Performance 


CONFIG. 

RUN 

PILOT 

COMMENT 

C-H 

RATING 

RMS 

PITCH ERROI 


5/7/1 

B 

- slight tendency to bobble 
in pitch and roll 

2.5 

0.89 

lC-11-1 

5/8/1 

B 

- could perform task easily 

- tried to be more aggressive, 
led to diminishing returns, 
possibly in worse longitudinal 
tracking scores 

2.5 

1.00 


5/31/1 

C 

- had to work a little 
-noticed pitch sensitivity 

3 

0.95 


5/7/3 

B 

- have to rest right arm on 
right leg for precise inputs 

- let oscillations die out, but 
could be aggressive 

5 

1.34 

1C-11-2 

5/29/5 

B 

- not putting smooth inputs in 

- really aggressive (noticably 
more aggressive inputs) 

- acquisition easy; fine 
tracking hard 

- hard to control this display 

- not as "tuned up" today, 
referring to piloting technique 

7-8 

1.33 

1C-11-3 

5/8/2 

B 

- applied learned technique, 
part of his compensation 

- flying aggressively but 
smoothly, did not let dot 
completely vibrate freely by 
"smoothing" 

- the performance was as good 
as previous run but considerable 
compensation required 

4 

1.13 
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pilot’s technique in the specified task. The comments associated with the run which received 
a rating of 8 imply that the pilot’s compensation was inappropriate, that is he was not “tuned 
up” properly. On the other hand, the comments for the run which received a rating of 4 imply 
that the pilot’s compensation, though complicated, was well chosen for the task. 

The above discussion indicates that the differences in the vehicle dynamics and in the 
experimental conditions can noticeably alter the properties of the inputs that the pilot 
produces. In addition, even though the RMS values for longitudinal stick inputs showed little 
variation due to the changes in the level of flexibility and the display dynamics that were used, 
the frequency characteristics of those inputs differ considerably. 


Simulator Dynamics 

To assure that the motion cues which are produced in such a simulation are improving 
the realism and not corrupting it, some knowledge of the simulator capabilities and 
experimental validity are needed. Frequency response data for the simulator was obtained to 
more completely document the simulation facility and address these points. This data was 
used, for example, to investigate the degree to which motion cues are attenuated by hardware 
limitations and signal filtering, thus affecting the validity of the simulation. Though 
simulation validity was not a specific objective of this study, the data that was collected could 
be used in subsequent analyses to address this issue. 

Figure 30 depicts a functional block diagram of the simulator. Frequency response data 
was collected at the locations identified by the letters A through E. The responses measured at 
these points are given in Table 20. The units and sign conventions for these responses are 
given in Table 21. The measurements were made for eight frequencies for the longitudinal 
responses (viz. to = 0.25, 0.5, 0.75, 1.0, 1.25, 1.5, 1.75, and 2.0 Hertz) and at three frequencies 
for the lateral responses (viz. co= 0.25, 0.5, and 0.75 Hertz). The amplitude of the force 
inputs was twenty pounds. 

The vehicle model that was used corresponds exactly to the model that was used in the 
simulation study. This includes the equations of motion, control system and control loader 
parameters. The frequencies of the two modeled structural modes were varied to address the 
effects of reduced stiffness on motion capability. Frequency response plots for the transfer 
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Figure 30 - Simulator Functional Block Diagram 









Tabic 20 - Frequency Response Inputs and Outputs 
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Table 21 - Scale Factors and Units 


Parameter 

Units 

Sign Conventions 

Fpitch* ^roll 

lbs 

Fpitch - positive aft 
F ro u - positive right 

^pitch’ ^roll 

inches 

0 - positive up 

0,p 

rad 

sec 

p - positive right 

q 

rad 

q - positive up 

2 

sec 

P 

rad 

sec 2 

p - positive right 

n z cp » "Zcaen 

g’s 

n Zq,’ "z**, ‘ P°si" ve down 

n ycp’ fiycoo, 

g’s 

"yep’ "yearn ' P° SitiVe ri g ht 

"ze*’ "yep 

g’s 

n^ - positive down 
n y<v - positive right 
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functions described in Table 20 are shown in Figures 31 through 35. These responses 
correspond to the baseline vehicle (i.e. co y = co z = 2.0 Hertz). Complete frequency response 
data for the vehicle model with several different vibration frequencies is presented in 
Appendix 6. 

As the structural stiffness is reduced, the responses associated with the vehicle dynamics 
(i.e. the responses between A or B and C, D or E) change. Note also that the data includes the 
time delays associated with numerical integration and conversion from digital to analog 
signals and vice versa. In addition, the washout filters and motion hardware are nonlinear, and 
so the frequency responses of the simulator are input amplitude dependent, in general. 

The result that the presence (or lack) of motion cues did not seem to significantly affect 
the simulated vehicle responses or the task performance, discussed previously was surprising. 
One would think that motion cues could be beneficial in providing additional information that 
the pilot can use in performing the required task. However, the pilot would also be inclined to 
limit his inputs in order to improve the ride quality and accept reduced performance when 
motion was present while becoming more aggressive and trying to improve performance when 
motion was not present. More aggressive behavior would result in larger inputs and lead to 
higher normal accelerations. But the RMS normal accelerations (shown in Figure 25) are not 
strongly dependent on motion status. 

This may have been due to several factors relating to the motion system of the simulator. 
Due to the limited travel of the motion base, washout is needed to avoid constantly hitting the 
displacement, rate or acceleration limits. As shown in Figure 40, the washout filters severly 
attenuate the motion cues, especially acceleration cues, thus reducing the amount of motion 
that the pilot feels. In addition, any bandwidth limitations of the motion system could make 
the accuracy of simulating the effects of structural vibration questionable, due to attenuation 
of higher-frequency motion and phase shifts which can cause the motion cues to be out of 
phase with the visual cues. This was in fact also indicated, for example, in Figure 40, 
especially in pitch acceleration. 
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Simulator Frequency Response 
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CHAPTER 4 

SUMMARY AND CONCLUSION 


The goal of this study was to develop an elastic aircraft model and conduct a simulation 
experiment to address the effects of aeroelastic dynamics on aircraft responses. This goal was 
accomplished with the vehicle model developed to describe the general motions of a 
representative large elastic aircraft and the simulation experiment documented in this report. 
This report included a discussion of all major aspects of the vehicle, including the equations of 
motion, aerodynamic data, geometry and structural data, control system dynamics, control 
loader data, and engine dynamics. In addition, all the information needed to implement the 
vehicle model on the Langley Visual/Motion Simulator was presented. 

A description of the experimental set-up and procedure, and a presentation of the 
experimental results were also included in this report. The task definitions, display dynamics 
and symbology, and configuration definitions were discussed. The results included pilot 
comments and subjective ratings, and important vehicle-response statistics and stick spectra. 
Frequency response data for the simulator were also included. 

A review and analysis of some of the pilot rating data and vehicle-response time- series 
data is also included. The results indicate that there are fundamental differences between the 
dynamics of rigid aircraft and elastic aircraft and that these differences can manifest 
themselves in increased tracking errors, degraded handling qualities, and changes in the 
frequency content of the pilot inputs. The results also confirm that handling qualities and 
tracking performance can significantly degrade as the level of flexibility is increased. The 
results also indicate the importance of the display dynamics. The definition of the displayed 
variables, or visual cues, have been shown to significantly affect results. 

Clearly, the reduction of structural stiffness combined with flexible display dynamics 
(those that incorporate the effects of structural deformation) lead to considerable degradations 
in handling qualities. The utilization of “rigid” display dynamics, which ideally filter out the 
effects of structural deformation, however, lead to little degradation in tracking performance 
regardless of the level of flexibility. That is, with the rigid display, results are less sensitive to 
variations in the level of structural stiffness than for the flexible display. This can be 
interpreted in terms of the pitch-rate-to-stick-deflection transfer functions listed in Chapter 2, 
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Table 7. For the rigid display, the effective pole- zero cancellation in the pitch-rate transfer 
function is noted, while in the flexible display case similar pole-zero cancellations do not 


occur. 
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CHAPTER 5 
RECOMMENDATIONS 


The general nature of this experiment and the variety of recorded response data that is 
available allows many issues which affect pilot/vehicle interactions for elastic aircraft to be 
addressed. The experiment was specifically designed to help study, evaluate, and develop 
predictive tools for pilot/vehicle performance, and to identify the degree and manner in which 
the simulator itself contributes to the pilot/vehicle interactions. The recorded time histories of 
all important pilot inputs and system responses are available for detailed statistical and time 
series analysis. In addition, the simulator model was developed to be readily altered and 
extended for use in future studies. Some specific applications of the experimental data and 
simulation model that could be performed in the future will be discussed. 

The response data from the pilot and vehicle simulation can be used to help extend and 
validate pilot modeling methods. For example, models for human operators need to be 
extended to account for multiple loop closures in multiple-axes tasks, and the experimental 
data is available to validate these extensions. The data can also be used to determine pilot 
describing functions that could then be used to identify the effects of motion cues and display 
dynamics on pilot dynamics. 

Frequency responses of the simulator were also collected. These can be used to further 
address the effects of simulator characteristics. For example, a model of the simulator could 
be developed and used in conjunction with the vehicle mathematical model and a pilot model 
(the OCM or crossover model) to quantify the effects of delays, lags, washout of motion cues, 
and other simulator dynamics on pilot/vehicle performance. Once these effects are identified, 
it may then be possible to better extrapolate results from experimental simulator studies to 
actual flight characteristics. In this way simulation could become a more valuable analysis 
and design tool. 
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The general, open structure of the elastic vehicle model that was used in the experiment 
has many applications for future research as well. Since it is a full multi-input, nonlinear, six 
degree of freedom model of a flexible aircraft and since it was developed to be easily altered 
and extended (for example, changing the level of stiffness, adding additional elastic degrees of 
freedom or implementing additional control systems) it can be used to study many aspects of 
advanced flight control development. Current flight control research is involved with 
multivariable control systems, robustness issues, functionally integrated control systems, etc. 
All of these areas could be studied experimentally using the elastic aircraft model developed 
for this study. 
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